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ABSTRACT 

Two  new  transient  testing  techniques  were  evaluated;  the  centrold 

method  developed  by  Kohlmayr  and  the  time  zero  intercept  technique.   The 

zero  intercept  method  was  found  to  be  the  most  promising  of  the  two  but 

is  limited  to  values  of  N   <  2.5.   The  centrold  technique  can  be  used 

effectively  when  the  value  of  N   is  less  than  5.0. 

tu 

A  heater  system  made  of  .001  inch  diameter  nichrome  wire  was  designed 
and  tested  to  determine  its  effect  on  the  transient  testing  of  matrix 
type  heat  exchangers.   Because  the  design  showed  no  improvement  in  the 
test  results  and  was  unreliable  its  use  was  discontinued. 
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NOMENCLATURE 


English  Letter  Symbols 


A 
A 


fr 


b 
C, 


H 


E 
G 

8c 
g 

h 


Kg) 

k 

k 


Matrix  tottil  heat  transfer  area 

Matrix  ii.nl.-iti  i  free  flow  area 

Matrix  total  frontal  area 

Matrix  solid  cross-sectional  area 
available  for  thermal  conduction 

Fin  thickness 

Flow  passage  perimeter  (  jS  Af  ) 

Fluid  stream  thermal  capacity  rate 
(mcf) 

Fluid  specific  heat 

Matrix  thermal  capacity  (W  c  ) 

s  s 

Fluid  specific  heat  at  constant 
pressure 

Matrix  material  specific  heat 

Flow  passage  hydraulic  diameter 
(4rh) 

Friction  power  per  unit  area 

Flow  stream  mass  velocity  (m/A  ) 

c 

Proportionality  factor  in  Newton's 
Second  Law 

Normalized  fluid  temperature  at 
inlet  to  test  section 

Surface  heat  transfer  coefficient 
for  convection;  heat  transfer 
power  per  unit  area  per  degree 
temperature  difference 

Deviation  from  step 

Fluid  thermal  conductivity 

Matrix  thermal  conductivity 


sq 

ft 

sq 

ft 

sq 

ft 

sq 

ft 

ft 

ft 

Btu/(hr 

deg  F) 

Btu/(lbm  deg  F) 
Btu/deg  F 
Btu/(lbm  deg  F) 

Btu/(lbm  deg  F) 
ft 

hp/sq  ft 

lbm/(hr  sq  ft) 

32.2  (lbm  ft)/(lbf 
sec*) 

dimensionless 
Btu/(hr  sq  ft  deg  F) 


dimensionless 

Btu/(hr  sq  ft  deg  F/ft) 

Btu/(hr  sq  ft  deg  F/ft) 


m 
P 
P 

q 

R 

rh 

t 
t* 

u 

V 
m 

w 

s 

w„ 


Total  matrix  flow  length 

Mass  flow  rate 

Pressure 

Matrix  porosity  (A  /A,.  ) 

c   fr 

Heat  transfer  rate 

Gas  constant  (53.35-air) 

Hydraulic  radius  (A  L/A) 

Temperature 

Normalized  temperature 

Flow  velocity 

Matrix  volume 

Matrix  mass 

Fluid  mass  in  matrix 

Distance  along  flow  passage 
from  the  matrix  inlet 


Reduced  length  (N   x) 

Li 


Greek  Letter  Symbols 


/ 


A 

e 

CEMTp 


r 


Compactness    (A/V   ) 
m 

Ratio  of  orifice  diameter  to  pipe 

diameter  (d  /d) 
o 

Difference  or  change  (time, 
temperature,  distance,  etc.) 


Time 

Free  time 
Fluid  viscosity 
Centroid  coordinate 
Density 


ft 

lbm/hr 

lbf/sq  ft 

dimensionless 

Btu/hr 

(ft  lbf)/(lbm  deg  R) 

ft 

deg  F 

dimensionless 

ft/sec 

cu  ft 

lbm 

lbm 

ft 

dimensionless 

sq  ft/cu  ft 
dimensionless 


sec,  hr 
dimensionless 
lbm/hr  ft 
dimensionless 
lbm/cu  ft 
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Subscripts 

atm  Local  atmosphere 

ave  Average 

f  Fluid  (gas,  air) 

i  Initial,  inlet 

m  Matrix,  mean 

0  At  orifice 

s  Solid  (Matrix  material),  static 

STD  Standard  (temperature  and  pressure) 

x  Local  conditions 

1  Inlet  conditions  (upstream  of  matrix  and  heaters) 

2  Inlet  conditions  at  matrix  entrance 

3  Exit  conditions  at  matrix  outlet 

Dimensionless  Groupings 

f  Fanning  friction  factor;  ratio  of  wall  shear 
stress  to  fluid  dynamic  head 


i 


2/3 
Colburn  j-factor  (N  N     ).   This  factor  plotted  vs, 

Reynolds  Number  defines  the  surface  heat  transfer 


characteristics . 


^        Longitudinal  heat  conduction  parameter  for  solid 

material  (k  A  /mLc.) 
s  s     I 

^Time    parameter    (hA0/W  c    ) 
s   s 
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1.    Introduction. 

The  transient  test  facility  at  the  Naval  Postgraduate  School (NPS) 
has  been  in  operation  for  several  years.   Two  techniques  have  been  used 
at  the  NPS  facility  for  the  determination  of  heat  transfer  data.   One  of 
these  is  the  maximum  slope  technique  developed  by  Locke  [13]  which  uti- 
lizes the  temperature-time  response  curve  of  the  fluid  leaving  the  matrix 
after  a  step  change  in  the  fluid  inlet  temperature.   Howard  [6]  extended 
this  technique  to  include  the  effects  of  longitudinal  conduction.   The 
other  technique  is  the  cyclic  technique  developed  by  Bell  and  Katz  [4], 
In  both  of  these  techniques  the  heat  transfer  parameter  being  sought  is 

N   ,  which  is  a  dimensionless  heat  transfer  parameter  equal  to  the  ratio 
tu 

of  the  convective  heat  transfer  rate  from  a  solid  to  the  heating  capacity 
rate  of  an  adjacent  fluid,  i.e., 

where:  ' 

Ji        -   unit  conductance  for  convection  heat 
transfer  (BTU/hr  sq  ft  deg  F) 

ri        »  total  heat  transfer  area  of  solid  (sq  ft) 

HI       =  mass  flow  rate  of  fluid  (lbm/hr) 

yCf     m   specific  heat  of  fluid  (Btu/lbm  deg  F) 

It  is  known  that  the  maximum  slope  technique  is  unreliable  at  values 

of  N   less  than  3.5.  due  to  the  large  errors  in  N   associated  with 
tu  6  tu 

errors  in  the  determination  of  maximum  slope ,   Furthermore,  it  has  been 

noted  experimentally  that  for  values  of  N   near  2.0,  the  temperature 

tu 

response  curve  of  a  fluid  displayed  its  maximum  slope  at  approximately 

time  zero  on  the  trace,  and  that  this  value  of  maximum  slope  was  not 

suited  for  determining  N 

6   tu 
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It  was  not  until  the  work  done  by  Kohlmayr  [9,  10],  which  gave  the 
exact  analytical  solution  to  the  single-blow  problem  based  on  Hausen's 
[5]  mathematical  model,  that  the  reason  for  these  errors  was  fully  under- 
stood.  Later  work  by  Kohlmayr  [11,  12]  demonstrated  ways  in  which  to 
handle  this  problem  as  well  as  the  development  of  a  new  technique  which 
might  be  used  in  place  of  the  maximum  slope  technique  for  values  of  N 

<5.   Kohlmayr  also  developed  a  means  to  handle  the  single-blow  problem 
when  other  than  a  step  change  in  the  inlet  fluid  temperature  was  made. 

The  purpose  of  this  thesis  was  first  to  determine  experimentally 
the  actual  inlet  temperature  response  of  the  NPS  facility  and  to  use  this 
known  temperature  response  in  the  manner  suggested  by  Kohlmayr.   Addi- 
tionally, a  new  heater  system  was  made  to  try  to  more  closely  approxi- 
mate a  step  change  in  the  fluid  inlet  temperature  so  that  the  physical 
process  and  mathematical  model  might  more  closely  resemble  one  another, 
thereby  improving  the  results  using  the  maximum  slope  technique. 

Another  technique,  the  "time  zero  intercept  technique,"  used  by 
Wheeler  [17]  was  investigated. 
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2.    Summary  of  Theory, 

A.  Background 

The  single-blow  transient  technique  which  is  used  to  determine 
heat  transfer  data  for  a  porous  solid  originally  used  for  its  mathe- 
matical model  Hausen's  [5]  partial  differential  equation  system.   The 
method  involved  comparing  the  recorded  exit  temperature  of  a  fluid  passing 
through  a  porous  solid  which  had  previously  undergone  a  step  change  in  its 
temperature ,  with  a  computed  response  curve  based  on  the  solution  to  Hausen's 
equations 

Locke  [13]  has  shown  that  there  exists  a  unique  relationship 

between  the  maximum  slope  of  the  response  curves  and  the  number  of  heat 

transfer  units,  N   .   However,  Hausen's  model  did  not  include  the  effects 
tu 

of  longitudinal  conduction,   Howard  [6],  by  the  use  of  a  finite  differ- 
ence technique  with  the  digital  computer,  included  the  effects  of  longi- 
tudinal conduction. 

It  should  be  noted  here  that  in  all  of  Kohlmayr's  work,  which 
will  be  discussed  later,  the  effects  of  longitudinal  conduction  are  not 
considered 

B,  Theory 

The  basic  assumptions  in  the  single-blow  problem  are: 

(1)  Properties  of  the  fluid  are  temperature  independent 

(2)  Fluid  flow  is  steady 

(3)  The  porous  solid  is  homogeneous 

(4)  The  thermal  conductivity  of  both  solid  and  fluid  is 
infinite  in  the  direction  perpendicular  to  flow 

(5)  Thermal  conductivity  of  the  solid  is  zero  in  the 
direction  of  flow 
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Hausen's  original  differential  equation  system  [5]  was  based  on 
the  energy  balance  between  a  fluid  passing  through  a  porous  solid  and  the 
solid.   The  system  of  differential  equations  which  resulted  from  this 
development  are  as  follows: 


<*t*  ts  -    ts 


<*> 


c 

in  which  a  "reduced  length"  variable, 

&      thjq,L        iu~  L 

and  a  dimensionless  time  parameter, 

K4  rhLx/\ws^s  i 

have  been  used. 

In  the  above  equations: 

X.S  -  temperature  of  the  solid  (deg  F) 
"C^  =  temperature  of  the  fluid  (deg  F) 
0       -   time  (hr) 

^C$   =  specific  heat  of  the  solid  (Btu/lbm  deg  F) 
]/\/s      -   mass  of  the  solid  (lbra) 

/X.        -   distance  along  flow  passage  measured  from  inlet  (ft) 
/_   =  total  length  of  solid  (ft) 
lA/^   =  mass  of  fluid  entrained  in  solid  (lbm) 
Kohlmayr  [11]  has  modified  Hausen's  original  equations  by  introducing  a 
new  dimensionless  time  variable  called  "free  time," 


m&  g  _   W*A*  ^  (1) 


WsJCs  Ws  As       L. 


or 
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(2) 


For  convenience  Kohlmayr  introduced  two  constants, 

OC  =  *"■*■  and    /3  =     W+X, 

which  are  fixed  for  any  particular  experiment.   For  most  practical  pur- 
poses,  /S  Z 0  and  JJ-  ^&CG    . 

Hausen"s  modified  equations  are  then  given  as 

JtHz,A),tSfcA)    _  #^,^  (4) 

where  t*  and  t*  are  the  temperatures  of  the  solid  and  fluid  respec- 
tively; which  have  arbitrarily  been  normalized.   The  above  equations  are 
subject  to  the  boundary  and  initial  conditions: 

*  z  -o  tr&jA.)  =  dC/L)  (5) 

and 

at/^-0        t?  (2,/L)       =    /  (6) 

where  g( yLL  )    represents  the  normalized,  time-dependent  fluid  tempera- 
ture at  the  inlet  cross-section.   Kohlmayr  solved  equations  (3)  through 
(6)  by  means  of  a  double  Laplace  Transform  [10]. 
The  results  of  Kohlmayr 's  solution  are: 


+ 'i  -  e  8  D  "fiUhX.  d^jy)  e^d  v] 


(7) 
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^~  .     is  an  entire  family  of  functions  of  the  order  K  (lC— 0, 1  .  .n^) 


CO 


00  = 


X 


niCn-v^l 


Evaluated  at 


n=o 


^■1      =  K 
equation  (7)  becomes 


-*L 


**(¥**)  r.hc  l-tfr* 


H- 


/C-/C^))e 


-a^CT-/*.) 


(7  -  <^(i/))oLi/ 


(8) 


(9) 


By  differentiation  of  equation  (9)  with  respect  to  Z/  , 


+  /£ 


(10) 


C  ^M(M-V)  -H,(Vt>v;)]  e^'^MdyJ 


With  these  results  Kohlmayr  [12]  then  showed  that  for  other 

than  a  step  change  in  inlet  temperature,  both  the  maximum  slope  and  "free 

time"  at  a  given  value  of  N   may  be  multi-valued.   He  further  demonstrated 

tu 

that  by  knowing  what  the  inlet  response  is,  one  can  determine  which  of 
the  values  of  maximum  slope  is  valid  and  also  how  to  determine  what  the 
error  bounds  are  in  the  use  of  some  known  inlet  temperature  change.   One 
of  Kohlmayr' s  results  is  that  the  correct  value  of  "relative"  maximum 
slope  occurring  in  the  fluid  exit  temperature  response  will  be  that  value 
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which  occurred  last  in  cases  in  which  KohlraayrJs  technique  is  applicable. 

Furthermore,  the  maximum  slope  method  is  unstable  for  N    <C   2.0,  is 

tu 

singular  at  N   =  2.0,  and  inaccurate  for  N   between  2.0  and  3.0.  [9,12] 
°         tu  tu  * 

In  view  of  the  short-comings  of  the  maximum  slope  technique, 
Kohlmayr  developed  an  indirect  curve -matching  technique  [11].   Some  im- 
portant results  of  this  development  are  summarized  below. 

With  the  solution  to  the  single-blow  problem  known,  it  is  pos- 
sible to  develop  an  indirect  curve-matching  technique  based  on  the  first 
moments  of  the  fluid  transient  response  curves.   This  method,  known  as 
the  centroid  method,  involves  the  reducing  of  both  the  theoretical  tempera' 

ture  response  t*(N   ,  JUL)    and  experimental  response  t*    (N   ,  M~)    into 
f   tu V fexp   tu  / 

two  different  single-valued  functions  based  on  the  one  parameter  N 

tu 

To  generalize  the  problem  further,  a  mapping  functional  was  defined, 
based  on  the  fluid  inlet  temperature  change  and  the  fluid's  exit  tempera- 
ture response  such  that 

where  the  following  restrictions  are  imposed  on  i^; 

(1)  ys     must  be  real,  single-valued,  continuous  with 
respect  to  both  t*(N   ,  M-   )  and  g( /C  ),  and  monotone  with  the  parameter 

tu 

(2)  J&   (N   )  must  be  monotone  increasing  with  1(g)  = 

Cy*~  s    the  "deviation  from  step." 


f?& 


(3) 


K.  >  where  K  is  some  measure 


cLAtt 

for  the  maximum  permissible  amplification  of  errors. 

(4)  For  any  N   and  any  given  deviation  from  step,  Y-   (N   ) 
tu  tu 

should  be  insensitive  with  respect  to  local  variations  of  g(X-<-). 


19 


(5)  The  evaluation  of  )£"   (N  )  must  be  simple  and 
straightforward. 

The  functional  chosen  was  the  first  moment  of  the  difference 
between  upstream  and  downstream  fluid  temperatures, 


,co 


</*  ^^/j  -^o^  =    /+    /     +/#*)<$*,  (id 


Z   '  M 

tcA. 


Due  to  the  difficulty  involved  with  integrating  equation  (11)  up  to 
values  of  JJ-   =  Oo    ,  a  new  functional  was  chosen  and  defined  as 


/>'°kv^^;-^j]°^  =  X+77-  +f}^d^     <i2> 


Some  further  restrictions  imposed  upon  g(  J^~)   are: 

(1)  g(  JJL   )  must  be  non-negative  monotone  decreasing 

(2)  have  initial  value  g(  O  )    =1  and 

(3)  assume  value  zero  for  all  free  times  which  exceed 

the  maximum  permissible  deviation  from  step:  %(M-   ) 

■  0  for  all  n      >     I 

yS~^      —       max 

Therefore  the  moment  functional  was  defined  in  terms  of  the  centroid  of 


the  area  under  this  difference  curve: 


A 


CEAJT£> 


—         so 


'%  (4*>j*)  -  tfri\<&L 


(13) 


V*& 


70 


£f(Miu_,/i)-p^)  <£/<_ 
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For  the  case  of  g(  JjL  )   ~   0  (step  change),  the  results  of  equation  (13) 
are  listed  as  Figure  1. 

In  order  to  use  the  above  results  for  any  given  inlet  tempera- 
ture change  g(  M    )t    Kohlmayr  arrived  at  the  following  empirical  relation- 
ship: 


A-P^**  *AJ*-*}  +  •" ^  ±  ■ 0I      (14) 


This  relationship  was  obtained  by  systematically  applying  the  centroid 
method  to  many  different  upstream  fluid  temperature  changes. 

Kohlmayr's  solution  to  the  single-blow  problem  may  also  be  used 
directly.   Returning  again  to  equations  (9)  and  (1) 


t:Wt,,ji)  =  /-  e 


-v.. 


>-n 


X) 


+('\2  2.,(fC(i/-j4  <?~'/:"  (1//L)(i-2  &))4  v 


(>. 


where 


u  =     WA  O 


W,  j£s 


")C 


t£(N  ,  M-   )  can  be  evaluated  at  "free  time,"  M-   =  0,  which  would  corres- 
pond to  the  time  that  an  element  of  fluid  which  has  undergone  a  change  in 
temperature  upstream  arrives  at  the  exit  of  the  solid.   Thus, 


' 


/  -  Q(o) 


f 


(15) 
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or 


/  -  tyA/,,0)  =  e^ 


[  i  -  9(0)] 


(16) 


Taking  the  natural  logarithm  of  both  sides  of  the  above  equation,  yields 


A> 


I  -  t*M  ,  o) 


l-QCo) 


v»- 


(17) 


Therefore 


-iW 


i-t;Hu,°) 


+  JU 


l-#(o 


(18) 


Note   that   for  a   step  change   in  the   inlet   temperature,    i.e.,    %{LL)   ~  0 
for  JLL>Ot    that  equation    (18)    reduces   to, 


tu 


i-t*(ti  ,o) 


(19) 


In  order  to  utilize  equation  (19)  directly  it  would  be  necessary  to  achieve 
a  step  change  in  the  inlet  fluid  temperature,  g(  LK  ).   This  would  occur 
at  real  time,  Q    =  0.   Then  measure  the  temperature  of  that  same  element 
of  fluid  which  had  undergone  the  step  change  in  temperature,  as  it  left 
the  porous  solid.   This  would  occur  at  "free  time,">U-  =  0. 
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Because  of  the  finite  response  time  of  the  different  components 
of  the  experimental  equipment,  i.e.,  heaters,  thermocouples  and  recorder, 
it  is  not  possible  to  meet  precisely  either  of  the  above  requirements. 
In  view  of  the  above  limitations  it  is  necessary  either  to  make  some  ap- 
proximations in  the  interpretation  of  the  temperature  response  data,  t*, 
or  else  devise  a  means  by  which  the  actual  response  data  might  be  treated 
directly.   Some  of  the  difficulties  involved  with  the  second  alternative 
are:   Refer  again  to  equation  (9) 


t?0L&  =    i-e 


MC 


i  -9C*> 


4 


A 


-c 


?/^VyU))  e^^'^f/-^))  cti/ 


which  may  be  stated  as 

+  J41  ^ML^-A))  e  (/-  $(jt) )  dv 

o 
and  the  expression  for  the  inlet  fluid  temperature, 

where  QH   is  the  equal  to  the  heater  time  constant  which  was  derived  in 
Appendix  C. 

One  approach  to  using  equation  (9)  would  be  to  measure  t*  at 
some  value  of  "free  time"  Ul  >   0  such  that  g(  LL    )  =  0  for  which  N 
would  be, 
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The  difficulty  here  is,  that  in  order  to  evaluate 
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N   must  first  be  known.   This  would  result  in  a  direct  curve  matching 
tu 

technique.   A  second  approach  would  be  to  measure  t*  at  some  time  //^>0 
such  that     /M- 

again  N   must  first  be  known.   In  both  the  above  methods  if   //C  is  small 
the  transient  response  of  the  sensing  and  recording  equipment  will  still 
present  problems. 

Because  of  the  difficulties  encountered  with  the  above  two 
methods,  the  first  alternative,  that  of  making  some  approximations  in  the 
interpretation  of  the  measured  response  data,  t*,  will  be  considered. 
Since  the  response  of  the  thermocouples  is  faster  than  that  of  the  .003" 
diameter  heaters,  and  approximately  equal  to  that  of  the  .001"  heaters, 
see  Appendix  C,  it  might  be  reasonable  to  assume  that  the  initial  response 
recorded  during  a  particular  run  is  that  response  due  to  the  heaters,  see 
Figure  2B.   Therefore,  it  is  possible  to  extrapolate  back  to  "free  time" 
yM-  -   0,  from  the  time  on  the  response  curve  where  the  transient  response 
of  the  heaters  has  died  out  and  treat  the  intersection  of  the  extrapolated 
curve  and  "free  time," JtA  =  0  as  an  actual  step  change  in  temperature.   The 
above  technique  is  referred  to  as  the  "zero  intercept"  technique. 
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3.    Experimental  Technique. 

The  existing  apparatus  at  the  NPS  facility,  see  Figures  3  and  4, 
has  been  designed  to  conform  to  the  idealizations  required  by  Howard  [5] 
for  the  use  of  his  conduction  parameter  in  the  maximum  slope  technique. 
Howard's  conduction  parameter^  ^  ,  is  defined  as 


where : 

-'Ks    »  thermal  conductivity  of  the  solid  (Btu/hr  sq  ft  deg  F/ft) 
r\^        ■  matrix  solid  cross-sectional  area  available  for  thermal 
conduction  (sq  ft) 
Howard's  idealizations  are  as  follows: 

(1)  The  fluid  flow  in  the  matrix  is  both  steady  and  uniform 
in  velocity  and  temperature  at  any  cross  section 

(2)  The  matrix  thermal  conductivity  is  finite  in  the  direction 
parallel  to  fluid  flow  and  infinite  in  the  direction  normal  to  flow 

(3)  The  matrix  thermal  conductivity  is  large  in  comparison  to 
that  of  the  contained  fluid 

(4)  The  thermal  properties  of  the  fluid  and  matrix  are  constant 
and  uniform 

(5)  The  convective  heat  transfer  coefficient  is  some  suitable 
average  and  remains  constant 

(6)  A  step  change  in  the  temperature  of  the  inlet  fluid  is  im- 
posed at  real  time  equal  to  zero. 

These  idealizations  result  in  making  the  fluid  flow  one  dimensional. 
They  in  no  way  conflict  with  the  restrictions  imposed  by  Kohlmayr  [11]. 

The  requirement  of  uniform  velocity  and  temperature  profiles  (1) 
is  met  by  a  specially  designed  entrance  nozzle,  flow  straightening 
screens,  and  an  even  distribution  of  heater  wires  across  the  channel. 
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Piersall  [15],  using  the  equipment,  verified  that  the  velocity  and  tempera- 
ture profiles  were,  indeed,  uniform  within  reasonable  limits.   The  im- 
portance of  uniform  velocity  and  temperature  profiles  was  demonstrated 
by  Wheeler  [17]. 

The  small  temperature  change  across  the  matrix  of  about  20°F  gives 
constant  thermal  properties  of  the  fluid  and  matrix  and  a  constant  heat 
transfer  coefficient  as  required  by  idealizations  (5)  and  (6). 

The  restriction  that  the  temperature  of  the  inlet  fluid  be  subjected 
to  a  step  change  will  be  discussed  later. 

The  data  required  for  the  computation  of  the  various  heat  transfer, 
fluid  flow  coefficients,  and  dimensionless  parameters  is  as  follows: 
•  atn\  =  atmospheric  pressure  (mm  Hg) 
y0  =  orifice  static  pressure  (inches  H-0) 

^p    ■  pressure  drop  across  orifice  (inches  H.O) 
A  p    =  pressure  drop  across  matrix  (inches  Ho0) 
Ps    ■  static  inlet  pressure  at  entrance  to  heaters 

(inches  HO) 
"t-o    ■  temperature  of  fluid  at  orifice  (millivolts) 
CLq         -   orifice  diameter 

/3    =  ratio  of  orifice  diameter  to  pipe  diameter 
CS    =  chart  speed  (sec/inch) 
U3-C/   =  downstream  response 
£3  "w.  ■  downstream  -  upstream  response 
Pressures  are  measured  with  either  a  draft  gage  or  a  water  manometer, 
depending  upon  the  orifice-flow  rate  combination,  except  for  atmospheric 
pressure  which  is  measured  by  a  mercury  barometer. 

The  temperature  responses  recorded  are  t«-t- ,  the  difference  between 
the  inlet  fluid  temperature  and  the  fluid  temperature  at  the  matrix  exit 
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(used  in  the  maximum  slope  and  zero  intercept  techniques)  and  t--t_, 
the  difference  between  the  fluid  exit  temperature  and  the  matrix  inlet 
fluid  temperature  (used  in  the  centroid  technique).   For  a  more  complete 
description  of  equipment,  the  reader  is  referred  to  Appendix  A.   See 
Figure  5  for  the  position  of  temperature  and  pressure  measurements. 

A  test  run  is  accomplished  by  predetermining  the  necessary  pressure 
drop  across  the  orifice  to  achieve  a  desired  flow  rate  and  to  determine 
the  number  of  heaters  necessary  to  achieve  a  20°F  temperature  rise.   Air 
is  drawn  through  the  apparatus  and  is  controlled  at  the  entrance  to  the 
turbocompressor.  When  the  desired  flow  rate  is  achieved  the  heaters  are 
then  energized  and  the  heated  air  and  test  core  are  allowed  to  reach  a 
steady  state  temperature,  at  which  time  pressure  measurements  are  re- 
corded  Power  to  the  heaters  is  then  secured  and  recordings  of  t~-t. 
and  t--t_,  both  as  functions  of  time  are  recorded  on  separate  channels 
of  the  Brush  recorder  (see  Figures  2A  and  B) .   The  temperature  at  the 
orifice  (t  )  is  measured  before  and  after  a  run  to  insure  the  same  ambient 
air  temperature  for  the  run. 

After  completion  of  the  desired  runs,  the  values  needed  to  compute 
slope  are  taken  from  the  recorded  traces  of  t.-t. .   The  maximum  value  of 
slope  is  obtained  visually  with  the  aid  of  a  straight  edge.   This  informa- 
tion is  included  on  the  data  sheet.   The  data  sheet  layout  conforms  to 
the  data  input  section  of  a  digital  computer  program  [18],  which  reduces 
the  data  used  in  the  maximum  slope  technique  and  calculates  the  parameter 

C?C   used  in  the  centroid  technique.  For  complete  details  on  data  re- 
duction the  reader  is  referred  to  Appendix  B  at  this  time. 
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4.    Description  of  Test  Matrices. 

The  two  cores  used  in  this  experiment  were  both  cores  which  had 
been  tested  previously  and  for  which  there  was  a  good  deal  of  data  avail- 
able. 

One  matrix,  the  Cercor  T20-38,  previously  tested  by  Howard  [7],  is 
a  ceramic  type-core  having  a  low  thermal  conductivity,  and  was  used  to 
eliminate  any  errors  that  might  be  introduced  by  longitudinal  conduction 
which  might  have  invalidated  results  using  the  centroid  and  zero  inter- 
cept techniques.   The  other  core  tested,  which  is  a  stainless  steel  plate- 
fin  type  matrix,  Solar  4,  had  been  tested  several  times  [15],  [16],  [18] 
and  was  used  so  that  comparisons  could  be  made  when  longitudinal  conduc- 
tion was  a  factor. 

Further  information  on  core  geometries  and  properties  is  shown  in 
Figures  6  and  7. 
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5,    Presentation  of  Results^ 

For  each  matrix  tested s  the  heat  transfer  and  flow  friction  character- 
istics have  been  computed.   The  computed  results  are  shown  tabulated  in 

Tables  I-II.   The  values  of  N   have  been  plotted  for  the  three  methods 

tu 

used  in  this  experiment.   Figures  8A  and  9  compare  values  of  N   based  on 

tu 

the  maximum  slope  and  centroidal  technique  for  cores  T20-38  and  Solar  4. 
Figures  8B  and  10  compare  values  of  N   based  on  the  maximum  slope  and  zero 
intercept  technique  for  the  same  two  cores.   Figure  11  is  a  comparison  of 
N   values  computed  by  the  maximum  slope  technique  for  the  two  heater 
systems.   Core  Solar  4  was  the  only  core  tested  with  the  new  .001  inch 
diameter  wire  heaters 

The  heater  response  curves  are  shown  in  Figure  12.   See  Appendix  C 
for  complete  details  used  in  determining  these  curves. 
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6.    Discussion  of  Results 

In  order  to  insure  that  the  technique  used  by  the  author  in  this 
experiment  was  correct,  the  results  of  the  single-blow  test  of  this  experi- 
ment were  compared  with  the  results  of  previous  tests  of  the  same  cores. 

Considerable  difficulty  was  experienced  in  the  use  of  the  new  (.001 
inch  diameter)  heater  system.   The  heaters  burned  out  frequently,  thereby 
limiting  the  range  of  Reynolds  Numbers  over  which  the  core  was  tested  and 
also  the  number  of  cores  actually  tested.   The  exact  reason  for  the  burn- 
out is  not  known.   It  was  initially  thought  to  be  due  to  surges  in  cur- 
rent caused  by  the  use  of  a  resistor  bank  in  the  first  experiments.  How- 
ever, when  the  resistor  bank  was  removed  and  voltage  to  the  heaters  control- 
led by  the  motor  generator  set  rheostat,  the  system  failed  again.   In  all, 
there  were  three  different  failures,  all  occurring  at  different  flow  rates 
and  with  a  different  number  of  heaters  in  use  at  the  time  of  failure. 

An  examination  of  Figure  11  shows  that  in  the  high  Reynolds  number 

ranges  there  is  little  difference  in  the  value  of  N   based  on  the  maximum 

tu 

slope  data  for  the  two  different  heater  systems.   Furthermore,  this  is  in 

a  range  where  the  errors  in  N   due  to  errors  in  the  maximum  slope  techni- 

tu 

que  are  greatest.  Without  actually  applying  Kohlmayr's  equations  in  the 

"extension  of  the  maximum  slope"  technique  [17],  it  is  impossible  to  get 

a  quantitative  value  for  the  error  in  N   due  to  deviations  from  the  step 

tu 

change.   However,  by  the  use  of  his  curves  based  on  the  "deviation  from 
step"  (previously  defined  as  1(g)  and  equal  to  the  area  under  the  inlet- 
temperature  response  curve)  which  are  presented  here  as  Figure  13,  it  can 
be  seen  that  for  I  <C  .  100  the  error  in  slope,  Am  which  is  the  difference 
in  the  maximum  slope  due  to  a  step  change  in  the  inlet  fluid  minus  the 
maximum  slope  due  to  a  non-step  change  in  the  inlet  fluid,  is  less  than 
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.01    It  must  be  pointed  out  that  Figure  13,  is  for  an  upstream  tempera- 

2 
ture  change  which  is  quadratic,  i.e.,  g(U   )  =  (1-M/A)  .   The  curves  of 

Am  versus  N    for  the  NPS  facility  would  be  slightly  different  from 
tu 

those  presented  in  Figure  13y  for  the  g(  LL)"s   associated  with  the  NPS 
facility's  heate:  systems  are  exponential.   However,  for  low  values  of 
I  (I  <^.  100)  the  error  in  maximum  slope  Am   would  be  small  regardless  of 
the  exact  form  of  the  inlet  temperature  response.   Since  the  maximum  1(g) 
resulting  from  either  heater  system  is  less  than  .03  it  would  not  be 
expected  for  one  to  detect  any  difference  of  slope  and  consequently  N 
resulting  from  the  use  of  either  system. 

Again  referring  to  Kohlmayr's  curves,  Figure  13,  deviations  from 
step  changes  in  temperatures  would  cause  large  errors  at  lower  flow  rates, 
i.e  ,  N   >5   However,  since  the  1(g)  of  the  NPS  facility  decreases 
with  decreasing  flow  rates  (see  Appendix  C)  it  could  be  assumed  that  the 
errors  in  maximum  slope  due  to  1(g)  would  be  small.   Furthermore,  since 
the  longitudinal  conduction  is  greater  at  low  flow  rates  its'  effects 
cannot  be  neglected  as  they  have  in  Kohlmayr's  assumptions. 

In  the  investigation  of  the  centroid  and  the  zero  intercept  techni- 
ques the  range  of  flow  rates  of  both  cores  tested  was  the  same  and  ranged 
from  nr\  ^   250  to  rn  ^   950  lbm/hr.   The  reason  for  950  lbm/hr  being 
the  maximum  value  tested  was  that  at  this  flow  rate  all  heaters  were  in 
use  and  in  order  to  go  to  higher  flow  rates,  an  Increase  in  voltage  would 
have  been  necessary  to  achieve  the  same  20  deg  F.  temperature  rise.   This 
would  ha  e  changed  the  heater  time  constant,  which  is  a  function  of  both 
flow  rate  and  the  number  of  heaters  in  use. 

Investigation  of  the  curves  of  Figure  8A  and  9  reveals  that  for  both 
cores  tested  there  is  very  close  agreement  between  the  results  predicted 
by  the  maximum  slope  technique,  shown  as  the  dotted  line,  and  the  results 
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of  the  centroid  technique  for  values  of  N   between  3.5  and  5.   In  fact 

n  tu 

there  is  even  some  overlapping  of  points  in  Figure  8A.   For  values  of 

N   below  3.5«  the  values  of  N   based  on  the  maximum  slope  technique 
tu  tu 

begin  to  decrease  from  the  predicted  values  whereas  the  values  determined 
using  the  centroid  technique  followed  closely  the  predicted  values  for 
the  case  of  Solar  4  and  increased  slightly  from  the  predicted  values  for 
core  T20-38. 

It  can  be  seen  from  Figure  8B8  which  compares  the  values  of  N 
computed  by  the  maximum  slope  and  zero  intercept  method,  for  core  T20-38, 
that  the  flow  rates  used  were  not  high  enough  for  the  zero  intercept 
method  to  be  used  to  its  best  advantage.   Only  those  runs  made  at  flow 
rates  greater  than  600  Ibm/hr  had  a  clearly  measurable  zero  intercept. 

Figure  10  again  compared  maximum  slope  N   's  with  those  evaluated 
by  the  zero  intercept  technique  in  this  case  for  core  Solar  4.   This  was 
done  for  both  sets  of  heaters.   It  is  important  to  note  the  close  correla- 
tion between  the  values  of  N   for  the  two  heater  systems.   This  will  be 

tu 

covered  more  fully  in  Section  7. 
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7.    Experimental  Uncertainties. 

Various  idealizations  and  boundary  conditions  have  been  imposed  for 
the  mathematical  model  of  the  physical  experiment.   Due  to  the  fact  that 
these  idealizations  and  boundary  conditions  have  not  been  precisely  met, 
certain  errors  have  been  introduced.   These  errors  are  difficult  to  as- 
sign a  numerical  value  to,  and  with  the  exception  of  the  deviation  from 
step  temperature  change  previously  mentioned  and  the  effects  of  longi- 
tudinal conduction,  will  not  be  discussed. 

The  experimental  errors  associated  with  the  maximum  slope  technique 
used  at  the  NPS  facility  have  been  covered  in  considerable  detail,  [2] [3] 
[15];  therefore  only  those  errors  introduced  in  this  particular  experi- 
ment will  be  discussed. 

In  using  Kohlmayr's  centroidal  technique  there  are  several  sources 
of  possible  error.   First,  Kohlmayr,  in  using  Hausen's  mathematical  model 
has  neglected  longitudinal  conduction  in  the  solid.   Howard  [6]  has  calcu- 
lated the  errors  in  N   associated  with  longitudinal  conduction  in  using 

tu 

the  maximum  slope  techniques,  but  no  such  information  is  available  for  use 
with  the  centroidal  technique,  for  the  analytical  solution  to  the  single- 
blow  problem,  in  this  case,  does  not  include  the  effects  of  longitudinal 
conduction.   Therefore,  one  can  only  get  a  qualitative  idea  as  to  the  ef- 
fects of  these  errors. 

From  Howard's  curves  of  N   versus  maximum  slope,  Figures  14  and  15, 
for  a  given  value  of  N   the  slope  of  a  temperature  response  curve  is 
different  for  different  values  of  the  conduction  parameter  A   .   There- 
fore, the  response  curves  and  consequently  the  centroid  of  the  response 
curves  would  be  different  for  different  values  of  A  .   The  difference  in 

maximum  slope  due  to   A  decreases  as  N   decreases.   For  example  at  N  = 

tu  tu 

30  a  difference  in  A.  of  approximately   .1  causes  a  907o  difference  in 
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maximum  slope ,  whereas  for  N   <  5  the  maximum  error  in  maximum  slope 

caused  by  a   .  1  difference  in  A  is  approximately  97o.   It  would  seem 

reasonable  to  assume  that  provided  Kohlmayr's  results  were  applied  only 

in  the  low  N   range  (N   <  5)  and  longitudinal  conduction  was  small  (i.e., 
tu        tu 

/\  *C  .1),  conduction  errors  could  be  neglected.   For  this  experiment 

the  maximum  value  of  A   for  either  core  at  values  of  N   <  5  was  less 

tu 

than  .022, 

Another  source  of  error  is  in  the  determination  of  the  centroid  of 
the  area  under  the  response  curve.   Two  methods  were  used  for  this  experi- 
ment, primarily  as  checks  on  each  other.   In  both  methods  the  critical 
point  is  the  determination  of  where  the  cut-off  is  for  a  particular 
curve.   The  cut-off  point  is  that  point  on  the  response  curve  which  is  one- 
tenth  of  the  maximum  value.   This  point,  defined  as  /^^'//O     occurs  at  a 
point  where  the  slope  of  the  curve  is  slowly  approaching  zero,  and  an 
error  of  one  millimeter  in  the  ordinate  can  cause  an  error  of  several 
millimeters  in  the  abscissa.   For  example,  in  the  data  taken  with  core 
Solar  4  when  N   =  1.41,  at  (tg-fe*)  =  6mm,  G  =  13.0  sec,  and  for  (tg-O 
=5,  C7  =13.75  sec.   This  was  for  a  typical  trace  in  which  (t^-t-)max  = 
50.5mm.  A  sketch  of  the  response  curve  (see  Figure  16)  will  help  to 
clarify  the  problem. 

The  actual  error  in  A^py^T^   caused  by  an  error  in  determining  i/Lu 
will  depend  on  the  value  of  N  .   A  typical  example  would  be  for  N   = 
3.5.  When  computed  by  the  manual  technique  (see  Appendix  B  for  details), 
it  was  found  that  a  +  1  mm  difference  in  ordinate  (which  is  &  2%   and 
would  be  a  maximum)  caused  an  error  of  3?  2%   in  /^pv^rn0      Ic  s*1011^  be 
pointed  out  that  these  values  are  approximate,  for  there  is  some  uncertain- 
ty associated  with  finding  the  intersection  of  the  three  centers  of  mass 
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lines,   shown  in  Figure  16  and  discussed  in  Appendix  B.   This  uncertainty 
usually  was  less  than  +  2%.   The  same  type  of  error  associated  with  find- 
ing the  cutoff  point  is  also  found  in  the  computer  technique  of  determining 
the  centroid,   In  this  method  the  computer  solution  required  data  points 
from  the  response  curve  at  fixed  intervals,  which  for  this  experiment  were 
5mm  intervals.   Therefore,  if  the  value  of  M-^/io occurred  at  a  point  mid- 
way in  any  one  5mm  increment,  there  could  be  as  much  as  a  2%  mm  error  in 
the  cut-off.   For  an  average  130  to  140  mm  trace  this  is  less  than  2%. 
For  comparative  purposes ,  data  runs  were  made  using  both  techniques  and 
the  maximum  deviation  in  centrold  between  the  two  methods  was  found  to  be 
4Z,  but  the  majority  of  deviations  was  less  than  1%. 

Another  difficulty  discovered  in  analyzing  the  temperature  trace  was 
caused  by  the  sensitivity  of  the  thermocouples.   During  a  particular  run 
a  change  in  the  steady  flow  caused  by  a  sudden  draft  in  the  laboratory 
caused  a  temperature  deviation  on  the  trace.   This  required  some  visual 
smoothing  on  the  part  of  the  author  to  extrapolate  the  actual  response. 
This  might  cause  errors  when  using  an  automated  data  reduction  process, 
unless  it  included  a  good  smoothing  technique. 

An  additional  problem  was  encountered  in  the  process  of  data  reduc- 
tion.  One  of  the  assumptions  made  in  using  the  centroid  technique  is  that 
at  time  zero  the  temperature  both  upstream  and  downstream  from  the  matrix 
is  constant.   It  was  found  that  this  was  not  the  case  for  this  experiment, 
the  reason  being  that  there  was  some  heating  of  the  thermocouples  located 
adjacent  to  the  heaters  due  to  radiation  from  the  heaters.   This  meant 
that  the  zero  reference  point  of  the  trace  of  (t»-t-)  established  before 
and  after  the  run  (heaters  off)  was  not  the  same  as  that  just  prior  to 
time  zero  when  the  heaters  were  deenergized  (see  Fig.  16).   Therefore, 
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where  the  response  curve  should  have  resembled  curve  a  of  Fig.  16,  it 

took  the  form  of  curve  b„   To  compensate  for  this8  point  c  was  taken  as 

time  zero  and  curve  a  used  to  determine  the  centroid.   It  should  be  pointed 

out  here  that  these  figures  are  greatly  exaggerated  here  for  clarity  and 

the  error  caused  by  this  difference  is  felt  to  be  negligible.  Also,  this 

difference  was  most  evident  at  low  flow  rates,  for  the  radiation  effect 

was  caused  primarily  by  the  heaters  nearest  the  thermocouples.   As  the  flow 

rates  increased  and  more  heaters  were  energized  the  percentage  of  heating 

due  to  radiation  became  less. 

So  far  no  mention  has  been  made  of  the  error  in  N   caused  by  an 

tu 

error  in  /A,^^.      Kohlmayr  [11]  has  conducted  a  linear  error  analysis 
y    CENTD 


based  on  the  approximation 

:  777s 


A  Hi 


and  introduced  a  relative  error  amplification  factor 

-ce/vr/) 


or 
A^ltuc      ~     Y   &/U- CENT& 

^i^  s^GE/JTD 

With  empirical   results  he  has   tabulated  values  of  K  versus     /^rv^Tn  anc* 

K  versus  N 

tu 

As  an  example,  using  Kohlmayr "s  curves  when  /^,_    is  measured  too 
high  by  21   and  when  K  =  -5.0,  then  N   as  obtained  by  this  method  would 
be  10%  too  low. 

Associated  with  the  zero  intercept  technique  are  two  distinct  types 
of  error.   One  results  from  the  initial  assumptions  on  which  this  technique 
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is  based.   These  assumptions  were  discussed  in  Section  2.   If  errors  are 
present  in  the  results ,  they  may  be  due  to  errors  in  the  original  assump- 
tions.  It  is  not  possible  to  get,  at  this  time,  a  quantitative  value  for 
this  type  of  error  or  in  fact  to  even  verify  that  the  technique  used  here 
is  valid.   The  only  factor  which  would  substantiate  the  technique's  vali- 
dity is  that  for  both  heater  systems  used  in  the  experiment  the  results 
in  terms  of  N   were  nearly  identical.   The  second  type  of  error  is  that 
associated  with  the  physical  data  reduction,  primarily  the  correct  extra- 
polation of  the  response  curve  back  to  time  zero,  see  Figure  17.   This 
may  be  seen  in  the  following  development: 

y^±  A/4,  =  -  At  0- ±;  ±    Ad/) 
e      e      *•        =  /- 1*  ±  a  t? 


7A^        _        j  ±      Ai* 


or 


±  a^     =  j^  a±  *t;  ) 

so  that  if  1-t*  =  .2  corresponding  to  a  value  of  N   =  1.609,  and 
f  tu 


or  27. 

or  1.25% 
or  -27. 


f-t? 

=  +  _ 

OZ 

th< 

»n 

**L= 

■*. 

02.02- 

if 

—  . 

OZ. 
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then 


A/1^=-.o/9g 


or   -1.27o. 
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8,    Conclusions, 

1.  The  present  design  of  the  new  .001  inch  diameter  heaters  was 
not  found  to  be  satisfactory  due  to  their  frequent  and  sometimes  un- 
explained fai lures ,   No  improvement  of  the  results  of  the  maximum  slope 
and  zero  intercept  tests  was  obtained  using  the  smaller  diameter  wire 
heaters.   No  results  were  obtained  with  the  new  heaters  for  evaluation 
with  the  centroid  technique ,  but  it  should  be  pointed  out  that  for  the 
centroid  technique  it  is  not  necessary  to  have  a  step  change  in  the  in- 
let temperature  so  long  as  the  actual  change  is  known. 

2.  It  would  be  extremely  difficult  to  try  and  compare  the  three 
techniques  used  in  this  experiment  on  the  basis  of  test  results  ac- 
curacy,  The  best  that  one  can  do  at  present  is  to  compare  the  various 
experimental  results  with  the  predicted  values  based  on  the  maximum 
slope  technique  evaluated  at  high  values  of  N    (N   >  5.0).   In  this 
respect  the  maximum  slope  technique  is  unreliable  for  N  <  3.5.   The 
centroid  technique  appears  to  have  good  results  for  N  <  5  to  as  low 
as  =75  which  was  the  minimum  value  tested  in  this  experiment.   For  both 
the  cores  tested  using  the  centroid  technique  the  results  were  either  the 
same  as  those  predicted  by  the  maximum  slope  technique  or  somewhat  higher. 
The  zero  intercept  method  appears  to  be  impractical  for  values  of  N   >2.5 
due  to  the  larger  differences  associated  with  the  logarithm  of  very  small 
numbers.   This  might  be  better  understood  by  looking  at  the  governing 
equation  in  the  zero  intercept  technique 


Mtu  =  -  jnO-t;^^) 


and  noting  that  as  the  argument  of  the  logarithmic  term  gets  smaller  the 
natural  logarithm  itself  approaches  negative  infinity.   Below  the  value  of 
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N   <2.5  the  results  were  higher  than  predicted  but  were  extremely  con- 
sistent, that  is,  for  different  runs  conducted  at  the  same  conditions 
the  results  were  nearly  identical. 

From  the  standpoint  of  ease  of  evaluating,  the  zero  intercept  techni« 
que  is  very  fast  and  simple  to  apply.   Furthermore,  in  the  method  used  in 
this  experiment  the  zero  intercept  technique  appears  to  be  insensitive 
to  slight  deviations  from  a  step  temperature  change.   On  the  other  hand 
the  results  of  the  centroid  technique  rely  entirely  upon  how  well  the  so 
called  "deviation  from  step"  change  of  the  fluid's  upstream  temperature 
is  known. 

The  centroid  technique  is  much  more  time  consuming  to  use  and  in 
view  of  the  fact  that  its  results  are  based  on  a  considerable  portion  of 
the  fluid  temperature  response  trace  it  is  more  subject  to  errors  caused 
by  sudden  fluctuations  in  the  ambient  conditions,  for  in  the  zero  inter- 
cept technique  only  the  very  first  part  of  the  trace  is  of  interest. 
Based  upon  the  above  conclusions  the  following  table  is  recommended  as 
a  guide  in  determining  which  of  the  aforementioned  techniques  should  be 
used  for  a  given  test  range. 


N   Range 
tu 

Technique 

3.5  <  N 

tu 

Maximum  Slope 

2.5  <  N  <3.5 
tu  ^ 

Centroid  Technique 

N   <2.5 
tu  ^ 

Zero  Intercept  Technique 
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9.    Recommendations  for  Further  Study. 

It  is  recommended  that  an  attempt  be  made  to  determine  both  experi- 
mentally and  analytically  the  effects  of  deviations  from  experimental 
assumptions  such  as  constant  fluid  properties  and  the  convective  heat 
transfer  coefficient  on  the  results  using  the  different  techniques 
including  the  cyclic  testing  technique. 
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APPENDIX  A 
Description  of  Equipment 
Air  Supply- 

The  working  fluid  for  the  experiment  is  air  drawn  through  the  test 
equipment  by  a  30  HP9  multi-stage  Spencer  Turbo-Compressor,  which  is 
rated  at  550  cfm  operating  on  a  220  V  a.c.  power  supply  see  Figure  3. 
Flow  Measuring  System, 

Flow  measurement  was  accomplished  with  an  ASME  standard  orifice 
section,,   Pressure  taps  were  located  d  and  d/2  diameters  upstream  and 
downstream  respectively  from  the  orifice.   Thin  concentric  orifices  with 
throat  diameters  of  .775,  1*232,  1.540,  and  2.310  inches  were  used. 
Heater  System. 

Two  different  heater  systems  were  used.   One  system  utilized  a  .0031 
inch  diameter  nichrome  wire  as  the  heating  element,  and  the  other  .001 
inch  nichrome  wire.   Both  systems  were  designed  to  give  the  same  nominal 
resistance  for  a  given  number  of  heater  switches. 

The  .0031  inch  system, (Figure  18A) ,  consisted  of  14  separate  bakelite 
frames,  each  wound  with  two  parallel-connected  heater  elements.   Each 
pair  of  heaters  was  controlled  by  an  individual  switch,  which  in  turn  was 
wired  in  parallel  with  the  other  heater  switches.   All  were  controlled 
by  one  master  switch. 

The  .0010  inch  system  (Figure  18B)  consisted  of  the  same  number  of 
bakelite  frames,  with  each  frame  containing  six  heater  elements  connected 
in  parallel.   The  number  of  heaters  and  voltage  could  be  varied  to  achieve 
approximately  a  20*F.  temperature  change  for  any  given  flow  rate. 
Matrix  Holder  and  Test  Section. 

The  matrix  holder  and  test  section  (Figure  4)  were  both  constructed 
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of  polyethylene  plastic.   The  matrix  holder  is  a  drawer  which  slides  into 
the  test  section.   A  removable  frame  on  the  side  of  the  holder  allows  for 
the  insertion  of  the  matrix.   The  flow  channel  is  3-1/16  inches  by  3-1/16 
inches,  and  can  hold  matrices  up  to  3  inches  long.   The  matrix  is  held 
snugly  in  place  by  styrofoam  insulation.   On  the  downstream  side  of  the 
holder  is  a  plate  containing  the  thermocouples  used  to  measure  t  .   The 
test  section  into  which  the  holder  slides  contains  the  heaters,  pressure 
taps,  and  the  thermocouple  set  downstream  from  the  heaters,  which  measures 

Inlet  Cone  and  Flow  Straightener 

This  section  was  designed  and  tested  by  Piersall  [15],  and  it  pro- 
vided a  uniform  velocity  profile  to  the  air  entering  the  matrix. 
Pressure  Measuring  System 

Pressure  taps  are  located  in  the  test  section  upstream  and  down- 
stream from  the  matrix  holder,  (see  Figure  5).   Two  other  taps  are  located 
at  the  orifice  section,   Each  pressure  tap  is  connected  by  flexible  tub- 
ing to  its  corresponding  manometer  and/or  draft  gage.   The  following  in- 
struments were  used  interchangeably,  depending  upon  the  flow  rate: 

1.  Ellison  Draft  Gage  Company,  0-3  in.  inclined  gage 

2.  Ellison  Draft  Gage  Company,  10  in.  manometer 

3.  Ellison  Draft  Gage  Company,  20  in.  manometer 

4.  Merriman  Instrument  Company,  120  in.  manometer 

5.  Precision  Thermometer  and  Instrument  Company,  mercury 
barometer 

Temperature  Measuring  System 

Temperatures  in  the  system  are  measured  at  four  locations:  (See 
Figure  5)inlet  to  the  system  (t..),  downstream  from  the  heaters  (t_) , 


45 


downstream  from  the  matrix  (t  ) 8  and  at  the  orifice  (t  ). 

3  ° 

Temperature  t„  is  measured  by  two  different  sets  of  30  gage  iron- 
constantan  thermocouples B   which  were  made  by  Traister  [16],   Each  set 
consists  of  5  thermocouples  connected  in  series.   Each  thermocouple  is 
individually  wrapped  in  teflon  tape  to  prevent  shorting.   All  ten  are 
contained  in  an  open-faced  aluminum  tube  mounted  in  a  frame  at  the  exit 
of  the  inlet  cone,   The  aluminum  tube  shielded  against  radiation  from  the 
heaters o 

Temperature  t  was  measured  by  a  set  of  five  „001  inch  diameter 
iron-cons tantan  thermocouples  connected  in  series.   The  output  of  these 
was  bucked  against  one  of  the  sets  measuring  t.  so  that  the  output  of  the 
two  sets  measured  (t.-t.). 

Temperature  t_  was  measured  in  the  same  manner  as  t_  and  the  output 

bucked  against  the  other  set  measuring  t.  so  that  the  output  measured 

1 

(t.j-t..).   The  outputs  (t_-t„)  and  (t.-t-)  could  be  recorded  separately  or 

bucked  against  each  other  to  give  (t--t_).   The  desired  outputs  were  then 

led  into  an  Astrodata  Model  886  Wideband  Differential  D.  C.  Amplifier, 

where  they  were  amplified  100sl„   From  there  the  signal  was  led  into  a 

Brush  Mark  280  Strip  Chart  Recorder. 

Temperature  t_  was  measured  by  a  30  gage  copper-cons tantan  thermo- 
u 

couple  referenced  to  an  ice  junction.   The  output  was  read  on  a  Leeds  & 
Northrup  Millivolt  Potentiometer. 
Heater  Power. 

The  power  to  the  heaters  was  supplied  from  the  250  V.  D.  C.  source 
in  the  laboratory.   The  actual  voltage  was  controlled  both  at  the  supply 
panel  and  at  the  apparatus  by  means  of  a  50  _fL  wound  resistor. 
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APPENDIX  B 
Data  Reduction  Relationships 

(1)  The  following  is  a  summary  of  the  data  reduction  relationships 
used  in  calculating  N   by  the  maximum  slope  technique. 
Geometry. 

Three  geometric  parameters  may  be  used  to  define  compact  heat  trans- 
fer surfaces.   This  allows  for  the  comparison  of  different  matrices. 

1.  Hydraulic  Diameter 

.     4x  free  flow  area  ,_  ,v 

H     h   heat  transfer  area 

2.  Porosity 

free  flow  area   =   c  (B-2) 

frontal  area       A, 

f  r 

3.  Area  Compactness 

S  ■  heat  transfer  area  =  A_  (B-3) 

'    matrix  volume        A,.  L 

fr 

Dividing  (B-2)  by  (B-3) 

rh  ■  p/|S  (B-4) 

Mass  Rate  of  Fluid  Flow 

The  mass  flow  rate,  m,  is  calculated  from  ASME  Power  Test  Code  [1] 
as  modified  by  Murdock  [14]  by  the  following  equation: 

m  =  359  Kc^fi  YjWr  (b-5) 


where 


c 

K  =        — —     flow  coefficient  including  velocity  of 
Y/  ffif  approach 

C  =  orofice  coefficient  of  discharge  [14] 

/3   =  ratio  of  orifice  diameter  to  pipe  diameter 

d  =  orifice  diameter  in  inches 
o 
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F  =  thermal  expansion  factor 

)f    "   P/RT  the  specific  weight  of  fluid  flowing,  assuming  a 
perfect  gas  (Ibf/fl 


Y  ■  expansion  factor 

P  =  absolute  static  pressure  at  orifice  (Ibf/sq  ft) 
R  -  gas  constant  for  air?  53»35  (ft-lbf/lbm°R) 
T  -  absolute  temperature  at  orifice  (deg  R) 
^P  =  pressure  drop  across  the  orifice  in  inches  H_0 
Substituting  the  expressions  for  K  and  ~£T  *n  equation  (B-5)  yields; 


no  =    3  5"5       c  do    £  V   lA£P 


•JTyS 


&T 


(B-6) 


From    [1]   Fig„    40A 


y  -      /    -  (oMl    -f  O.  35/6    )   X_ 

1  * 


o     jCf>  / 


xt 


/v- 


k  =  lo4  for  air9  rati 

Also,  from  [1],  Fig.  38 

F  =  1.0 
a 

P  =  (P   -P  /13o6){0o4912*144)  lbf/ft2 
atm  o 

P  .   =  local  atmospheric  pressure  in  inches  He 
atm  ° 

P  =  static  pressure  upstream  of  the  orifice  plate  in  inches  Hg 
o 

Making  the  above  substitution  in  (B-6)  with  the  constants  necessary  to  be 
dimensionaiiy  consistent  yields; 


m  =  539. 8/  (e/JT^JF)  dj 


l-(M 


-+  .035/3^) 


&Po 


(fa!)  n*i 


t0  t¥59j 


V/2. 


(B-7) 
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Reynolds  Number 

Reynolds  Number  is  defined  as; 

where  G.  the  mass  flow  velocity  ■  m/A  ■  m/pA,. 

c       fr 

and  U.     is  the  fluid  viscosity. 
Substituting: 

Mec  =  w  0>H  /Jul  fl^-P  =  ^rin/jjL  f)/r  B 


(B-8) 


(B-9) 


or 


//^     =  i+mL  / fijjL 


(B-10) 


Maximum  Slope 


The  maximum  slope  of  the  downstream  cooling  curve  (t_-t.)  is  a  unique 
function  of  N    [13]  and   At6!- 


tu 


t,  -t- 


O-  ,  N   and  ^  have  been  previously  defined,  but  for  convenience  are 


restated 


rruc* 


inh/:fL 


the  new  temperature  introduced  is  t.  and  is  equal  to  the  temperature 
of  both  the  solid  and  the  fluid  prior  to  cooling. 
Therefore: 


t/AL  =  0»4  /k/s  t, )  O 
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and 


letting; 


and 


cL  (2-/^ )  =  (nub/  t/s  a  )  do 


C..  -  inc.  (Fluid  stream  thermal  capacity  rate) 


C  -  W  c   (Matrix  heat  c 
s    s  s 

equation  (B-ll)  becomes „ 

cS(}/*4J  =  (o/c)  do 

Furthermore 

*'  -*i  tt  -±: 

whose  derivative  is; 


(B-ll) 


\t,  -±i 


I      oL  (i*  -t, ) 


t  Y-ti 


(B-12) 


combining  equations  (B-ll)  and  (B-12)  yields 


U,  -a 


KA       /      U(U-t,) 


<*ft/A4J   Lx    "^    t,-*i         do 


(B-13) 


At  =  t<  -t, 


(ti-t) 


-£irn<P    0 
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From  the  above  sketch  of  the  downstream  temperature  response  curve 


*fr3-t,) 


do 

X/  chart  speed 


Mfh 


d(t3-t1) 


combining  with  the  ratio  of, 


matrix  capacity 


flow  stream  capacity  rate 
and  equation  (B-13) 


=     Y/X 
=  d6 

=  Y 

-  A  t 


=  _s      (1/sec) 


cl, 


-t.-t: 


*(*K 


-      Cs        I       Y     Ch(Lri    Speed, 
^f    At     X 


(B-14) 


_  W*x 


This  value  of  maximum  slope  and  /[      are  then  used  to  enter  Table  III 
or  Figure  14  or  15  to  get  the  corresponding  value  of  N 
(2)  Centroidal  Technique 

For  the  centroidal  technique  the  recording  trace  of  (t.-t.)  is 
used.   Two  methods  are  actually  involved  in  utilizing  this  trace.   One 
method  is  to  copy  the  trace  physically  with  carbon  paper  on  to  a  thin 
piece  of  cardboard,   Then  the  trace  copy  is  cut  out,  and  by  means  of 
a  "plumb-line"  attached  to  a  pivot8  which  allows  the  trace  to  swing  free, 
the  line  passing  through  the  pivot  point  and  the  center  of  mass  is  determined. 
This  is  done  for  several  points  (at  least  3)  and  the  intersection  of  the 
different  lines  is  the  centroid  of  the  area  under  the  response  curve. 
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These  lines  are  shown  in  Figure  16.   Ideally  all  lines  will  intersect  at 

the  same  place 9  however 9  this  is  not  the  case  in  actuality  due  to  errors 

in  manipulation. 

An  alternate  method  is  also  used  to  find  the  centroid  of  the  function 

(t--t-)o   This  involves  the  use  of  a  computer  program9  which  when  fed 

actual  curve  data  (normalized)  computes  the  centroid.   This  is  done  in 

the  absence  of  an  apparatus  which  could  convert  directly  the  temperature 

response  to  information  useable  by  the  digital  computer,  such  as  a  paper 

tape  puncher. 

Once  M.„m-mn.   is  determined  from  the  trace  of  (t^-t-)  it  is  necessary 
/  CISNTB  J   2 

to  compute  1(g)  in  the  manner  described  in  Appendix  C.   This  value  of  1(g) 
is  then  used  to  determine  JLL0     .   Recall  Equation  (14) 


or 


Once     M~o     is  determined 9   Kohlmayr's  curve   of  N       vs.>U0  ,   Figure   1, 


may  then  be  entered  to  give  N 

tu 


Zero  Intercept  Technique. 
To  interpret  the  physical  data9  i.e.9  the  recorded  trace  of 
(t~-t.)9  which  represents  tf9  in  the  manner  described  in  Section  3,  of 
this  reports 

(i)  Extrapolate  by  means  of  a  French  curve  the  response 
curve  (t~-t_)  from  some  position  on  the  trace  where  the  transient  response 
of  the  heaters  has  died  out9  back  to  the  vertical  line  passing  through  the 
pointy^-  =  0.   This  point  would  be  the  position  where  the  first  change  in 
temperature  was  detected  by  the  thermocouples  and  is  referred  to  as  (t~-t  )0 
or  tf  (N£u  ,0  )9   Figure  2B.. 

(ii)  Record  the  value  of  (t«-t.)  at  the  intersection  of 
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the  time  zero  line  and  the  extrapolated  response  curve,  also  record  the 

value  of  (t_-t.)  at  the  time  just  prior  to  the  heaters  being  deenergerized. 

This  temperature  is  referred  to  as  (t0-t.) 

J   1  max 

(iii)  N    is  then  equal  to, 
tu 

Ntu  =  -too-  cvhV^VhW 

which  was  developed  in  Section  3,  as, 

A 


tLr-**.0-tf(*^fl)) 
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APPENDIX  C 

To  determine  the  transient  response  of  the  NPS  facility  for  differ- 
ent heater  sets,  runs  were  conducted  in  which  the  temperature  difference 
(t_-t-)  was  recorded  as  a  function  of  time.   This  temperature  difference, 
which  was  previously  called  "upstream  temperature"  is  the  difference  in 
temperature  of  the  air  entering  the  test  rig  and  the  air  leaving  the 
heaters. 

Analytical  determination  of  the  heater  time  constant  is  straight- 
forward and  presents  no  difficulty  based  strictly  on  the  analysis  of  a 
cylinder  in  cross  flow.   The  analytical  determination  of  the  thermocouple 
time  constant  is  more  difficult  and  requires  some  judgment.   Since  the 
thermocouples  are  constructed  of  .001  inch  diameter  wire,  it  is  not  pos- 
sible to  describe  precisely  and  mathematically  the  geometry  of  their 
junction.   In  the  construction  of  the  thermocouples  the  two  different 
wires  used  crossed^  at  an  angle  approaching  180  degrees,  and  arc-welded 
at  the  junction.   Therefore ,  the  junction  was  neither  spherical  nor 
cylindrical  (which  is  the  case  when  the  wires  are  butt-welded).   For  this 
reason  the  model  of  the  junction  was  arbitrarily  chosen  to  be  spherical 
and  calculations  were  based  on  a  spherical  junction  .002  inches  in  dia- 
meter.  On  this  basis  it  was  determined  that  the  time  constant  of  the  .003 
inch  heaters  was  nearly  three  times  greater  than  the  time  constant  of  the 
thermocouples 9  while  the  time  constant  of  the  .001  inch  heaters  was  ap- 
proximately two- thirds  that  of  the  thermocouples. 

Since  the  response  of  the  heater-thermocouple  circuit  would  be  the 
sum  of  at  least  two  exponential  terms 9  which  would  be  difficult  to 
separate  experimentally,  the  terms  have  been  treated  here  as  one  expon- 
ential term,  which  is  a  good  approximation  of  the  actual  system  response 
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if  one  neglects  the  initial  points  on  the  response  curve.   From  the 
normalized  temperature  response  curves  the  following  empirical  relations 
were  developed? 


t*  -t,      _  e 


-©. 


ft.-t> 


M*\* 


where  6/  =  time  in  seconds  and 

t/fy      is  the  heater  time  constant 
in  seconds  (C-l) 


k 


°3 


G      loc.  a 


Q, 


H 


(C-2) 


Oh  = 


-e 


Mf\\ 


,o;e 


(C-3) 


where  the  value  of 


-  e 


log 


(t  -t^max. 


is  the  reciprocal  of 


■  tie  5. ope  of  the  temperature  response  curve. 

The  heater  time  constant  was  then  plotted  against  mass  flow  rate  (  /^  ) 
From  these  curves  the  following  relationships  were  obtained: 


log  0/t     =  -  a  log  hr,     +  log  £ 


where  a  «£•  c  are 
constants 


(C-4) 
(C-5) 


Since  the  centroid  technique  makes  use  of  the  deviation  from  step  defined 


55 


1(9)    =      J  dr^)dM-  (c-6) 

o 

JUL     =   <*B  and       CtjUu^  <*de 

en         ^^-)    =   ^^     =     <?~     ^     =      ^f^^rf  (C-7) 

Z(9)=/  e  dp  =    *Q»  (c-8) 


as 

and   since 
th 


r<? 


where  again 


c/v  "  * =    rrl*       constant    (for  any  given 

core) 


Therefore  to  determine  1(g)  for  a  particular  run,  enter  the  curve  of  Q^ 
vs.  W)       (Figure  12)  with  appropriate  jfY\     .   Note  that  this  curve  is 
independent  of  the  matrix  used  in  the  run.   Then  multiply  £>/   by  the 
appropriate   ©C   to  obtain  1(g)  used  in  centroid  technique.   Figure  12 
shows  the  experimental  values  of  6>y/  vs.  DO        plus  the  theoretical 
plots  of  heater  and  thermocouple  transient  response  for  both  sets  of 
heaters. 
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Figure  9.    Core  Solar  J+  Teat  Results 
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